Melting of Moving Vortex Lattices in Systems With Periodic Pinning 
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We study numerically the effects of temperature on moving vortex lattices interacting with periodic 
pinning arrays. For low temperatures the vortex lattice flows in channels, forming a hexatic structure 
with long range transverse and longitudinal ordering. At higher temperatures, a transition to a 
smectic state occurs where vortices wander between channels and longitudinal order is lost while 
transverse order remains. At the highest temperatures the vortex lattice melts into an isotropic 
liquid. 

PACS numbers: 74.60.Ge, 74.60.Jg 



Driven vortex lattices interacting with quenched and 
thermal disorder are an ideal system in which to study 
the nonequilibrium phases and transitions that arise from 
the interplay of competing interactions ■ Transport 
measurements || , neutron scattering || , and Bitter dec- 
oration experiments Q have provided strong evidence 
for transitions between different vortex dynamic phases, 
including creep, plastic flow, and ordered (elastic) flow. 
Theoretical work and simulations Jlj] suggested that for a 
moving vortex lattice, at high velocities the effect of dis- 
order can be represented via a shaking temperature T s h, 
inversely proportional to the velocity. At high velocities 
T s h decreases below the melting temperature T m , and the 
vortices reorder. T m of the ordered moving vortex lattice 
is near the equilibrium T m of the disorder-free station- 
ary lattice. In the highly driven, solidified state, theo- 
retical H work suggested that the vortex lattice forms 
a moving Bragg-glass (MBG), or a strongly anisotropic 
moving smectic (MS) or moving transverse glass (MTG) 
H . In the MBG the vortices move in correlated channels 
with few defects, producing quasi-long range order. In 
the MS/MTG vortices move in uncorrelated channels, so 
although power law transverse order is present, the lon- 
gitudinal order is short range only. Strongly anisotropic 
ordering, consistent with a MS/MTG, as well as more 
ordered phases, consistent with a MBG and vortex chan- 
neling, have been seen in simulations and Bitter dec- 
oration experiments . 

Despite the considerable work done in the case of ran- 
dom quenched disorder, the effects of thermal disorder 
and melting for the interesting case of a moving vortex 
lattice interacting with a periodic pinning substrate have 
been far less studied. Periodic pinning substrates in su- 
perconductors can be created with arrays of microholcs 
H and magnetic dots j|. In all such systems, the pin 
radius is smaller than the distance between pins, so that 
a moving vortex spends the largest fraction of its time in 
the unpinned area. A crucial difference from the random 
pinning case is that the effect of the periodic pinning 
cannot be represented by a shaking temperature. The 
same periodicity also induces true long range transverse 
or longitudinal order in the moving lattice. 



We report a numerical study of the melting of mov- 
ing vortices in square periodic pinning arrays. With no 
driving, the system with pinning melts at a higher tem- 
perature than the pin free system. For moving vortices 
at low temperatures, we observe a triangular lattice flow- 
ing in strict ID correlated channels, with transverse and 
longitudinal long range order. The transverse order is 
greater than the longitudinal order, and the anisotropy 
increases with temperature. At higher temperatures near 
the melting temperature of the clean equilibrium system, 
we observe a transition to a moving smectic (MS) state, 
where transverse vortex wandering between channels oc- 
curs. At even higher temperatures near the melting tem- 
perature of the pinned equilibrium system, the moving 
smectic melts into a moving isotropic liquid (ML). We 
present the dynamic phase diagram and explain its fea- 
tures in terms of the pinned and unpinned equilibium 
melting transitions. 

We use finite temperature overdamped molecular dy- 
namics simulations in two dimensions, = ryv^ = 
f™ + fj Up + frf + fj T . We impose periodic boundary condi- 
tion in x and y. The force between vortices at and iv, 
is f.™ = J2fZi foA^Qvi ~ r^l/Ajfy, where K^r/X) is 
a modified Bessel function, /o = $q/87t 2 A 3 , A is the pen- 
etration depth and the parameter A v tunes the vortex- 
vortex interaction strength. For most of the results here 
A v = 1.0. The driving force f^, representing a Lorentz 
force, is in the x direction. The pinning is modeled as 
attractive parabolic traps of maximum strength f p and 
a range r p which is less than the distance between pins 
a. The thermal force ff has the properties < ff >= 
and < fi{t)fj{i!) >= 2r 1 k B T5 lj 5{t - t). The temper- 
ature T = (l/2rjkB)(Afo) 2 St where St is the time step 
in the simulation and A is the number we tune to vary 
T, with fj = Af - We take /o = &s = 77 = 1 and use 
St = 0.04. We explore the phase diagram by conducting 
constant fd or T sweeps on the fd — T plane. Our model 
is most relevant to superconductors with perodic arrays 
of columnar defects or thin-film superdonctors where the 
vortices can be approximated as 2D objects. A realistic 
3D model would be needed to address the exact nature 
of the liquid phase, such as whether it is a line liquid. In 
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FIG. 1. The structure factor peaks < S(k) > /N v versus 
T for a system with N v — N p . Upper curve: Sample with 
a square pinning substrate of strength f p /fo = 0.22. Lower 
curve: Sample containing no pinning. The melting for the 
pinned system occurs at much higher T than melting in the 
clean system. Inset: V x versus /<j with curves for f p /fo = 0.22 
and T = 0.0001 and 0.0725. 

this work we consider only the commensurate case where 
the number of vortices N v equals the number of pin- 
ning sites N p . Results for the incommensurate cases will 
be presented elsewhere. The initial vortex positions are 
generated by simulated annealing with each pinning site 
capturing one vortex. For finite size scaling we analyze 
system sizes with N v ~ L 2 for N v = 224 to N v = 2112. 

We first establish the melting temperature in the pin 
free and pinned systems without external drive (fd = 0) 
with the vortex displacements d r =< \(r(t) — r(0))| 2 >, 
and the structure factor S(k) = J2i, 3 ■ e^'^'^^l 
In the pinned system the vortex lattice has the same 
square symmetry as the pinning lattice. The melting 
temperature is determined from the simultaneous onset 
of diffusion and a drop in the peaks in S(k). In Fig. 1 we 
show that the melting temperature T m is higher in the 
pinned system, with TJ^ rs 0.03, than in the unpinned 
system, T£f » 0.0058. This is reasonable since at com- 
mensuration the pins stiffen the vortex lattice. 

Next we explore the dynamic phases of the system. For 
f p = 0.22/o the T — depinning occurs at fd = 0.22/ - 
For fixed fd — 0.45 we perform a T sweep, and monitor 
S'(k). Every 400 MD steps we measure the transverse 
displacements d y —< \y(T) — y(0)| 2 > from the initial 
positions of the vortices at T = 0. For low drives the 
vortices form a pinned square vortex lattice. The pinned 
phase is defined by measuring V x = (1/N V ) J^i=i v i 'X. In 
the inset of Fig. 1 we show the typical V x versus fd curves 
for two different temperatures. The transition from the 
pinned to moving phases are marked by a jump in the V x 
at a well defined fd- From V x versus fd curves we found 
little evidence for plastic or collective creep behaviour 
for temperatures below T^; however, much longer time 
scales would be necessary to explore the creep behaviour. 

Above some T dependent driving force, we find that 
the vortices form a triangular lattice, with the principle 
lattice vector aligned with the direction of motion, since 
the vortices spend part of thier time between pins where 
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FIG. 2. (a) The transverse structure factor peak 

S{kr)/N v (upper curve) and the longitudinal peak S(k]f)/N v 
(lower curve) for a 48A x 48A system with fd = 0.45/o, 
f p = 0.22 f p , and N v = N p . (b) The average transverse 
displacements of the vortices from their initial positions, 
d y =< \y(T) - 2/(0) | 2 >. Inset of (a): the scaling of S(k)/N v 
for different system sizes where N„ ~ L 2 . For T = 0.0025 
both the transverse (circles) and longitudinal peaks (squares 
and shifted down by 0.3 for clarity) scale as S(k) ~ L°; 
for T = 0.015 the transverse peaks (up triangles) scale as 
S(kr) ~ L while the longitudinal peaks (left triangles) scale 
as S(kL) ~ I/ -1 ' 96 . Inset of (b): the evolution of the fraction 
of six-fold coordinated vortices Pe versus temperature. 

the vortex-vortex interaction dominates. 

In Fig. 2(a) we plot the transverse S(kT) and lon- 
gitudinal S(fci) structure function peaks for increasing 
T. For low T, S(k T ) is only slightly larger than S(k L ). 
This anisotropy between the peaks becomes more pro- 
nounced as T is increased. The ordering can be ana- 
lyzed from the finite size scaling of the structure factor 
as S(k)/L 2 ~ L~'' |ll]]. A solid with long range order 
will have 77 = while a system with short range order 
such as a liquid will have 77 = 2. In the inset of Fig. 2(a), 
S(k T )/L 2 and S(k L )/L 2 (where N v ~ L 2 ) are plotted 
for different system sizes in the low temperature regime 
(T = 0.0025). Both peaks scale as 77 rs 0, indicative of 
long range order. Near T — 0.0095, which we label Tms, 
S(k[ J ) drops precipitously, indicating complete loss of lon- 
gitudinal order, while S(kr) retains a significant finite 
value. Gradually S(kx) drops to the level of S(]cl) near 
T = 0.03, which we label T M l- In this T > T MS regime 
the longitudinal peak S(kr,)/L 2 scales as 77 = 1.95 ±0.03, 
consistent with an r\ = 2 scaling behavior indicating the 
loss of longitudinal order. At this same temperature the 
transverse peak S(kT)/L 2 shows a ~ L~°-° form (trian- 
gles in upper inset of Fig. 2), indicating that long range 
transverse order is still present. The behavior of the two 
peaks for Tms < T < Tml indicates the presence of a 
moving smectic phase. For T > Tml both S{kif) and 
S(kx) scale as ~ L~ 2 as the system becomes an isotropic 
liquid. We label these three phases the moving crystal 
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FIG. 3. The real space vortex images (a,c,e) and vortex 
trajectories (b,d,f) for the same system as in Fig. 2 for the 
moving lattice (T = 0.0025) (a,b); the moving smectic (T = 
0.015) (c,d); and the moving liquid (T = 0.035) (e,f). 

'(MC), the moving smectic (MS) and moving liquid phase 
(ML). 7] has some temperature dependence near the tran- 
sitions which we will examine elsewhere. In the inset to 
Fig. 2(b) we plot the fraction of 6- fold coordination num- 
ber P 6 versus T as obtained from the Voronoi construc- 
tion. The proliferation of defects occurs at T = Tms as 
seen by the drop in P 6 . In Fig. 2(b) we plot the trans- 
verse displacement d v . For T < Tms, d y w indicating 
that the vortices are moving in straight ID channels. For 
T > Tms, d y increases indicating that the onset of vortex 
wondering in the y direction is correlated with the loss 
of longitudinal order and the proliferation of defects. 

To examine the individual vortex behavior in the differ- 
ent dynamic phases in Fig. 3 we plot snapshots of the vor- 
tex positions and trajectories for MC at T = 0.0025, MS 
at T = 0.015, and ML at T = 0.035. In the MC phase, 
the vortices form an ordered triangular lattice [Fig. 3(a)], 
and move in correlated ID paths [Fig. 3(b)], which agrees 



FIG. 4. The dynamic phase diagram of fd versus T. For 
open circles and filled squares, the current is fixed and the 
temperature is increased. For the filled circles, the tempera- 
ture is fixed and the drive is increased. The pinned to moving 
phases are determined from fd versus vortex velocity V x . The 
boundaries of the other phases are determined from trans- 
verse and longitudinal < S(k) > and d y . The open circles 
correspond to the MC to ML transition and the filled squares 
correspond to the MS to ML transition for a system with the 
same parameters in Fig. 2. The filled circles correspond to 
the transition out of the pinned phase P. Inset: The phase di- 
agram for A v , the vortex- vortex interaction prefactor, versus 
T. The triangles correspond to the zero drive clean melting 
temperature while the open circles and squares correspond 
to the MC-ML transition and the MS-ML transition with 
f d = 0.45/q. 

with the zero transverse wandering d y s=s in Fig. 2(b). 
As T is increased in the MC phase the channel width 
increases; however, vortices do not cross from one channel 
to another. The pinning also induces the long-range order 
seen from the scaling of S(k), whereas algebraic decay 
would be expected since the system is in 2D. 

In the MS phase, shown in Fig. 3(c), the vortices are 
much less ordered than in the MC phase; however, some 
order remains. The vortex trajectories in the MS phase, 
Fig. 3(d), reveal that although some channeling occurs 
along the pinning, there is considerable vortex motion 
between and across the channels. This interchannel vor- 
tex motion accounts for the increase of the transverse dis- 
placements, d y , at the onset of the MS phase in Fig. 2(b). 
The residual vortex channeling accounts for the finite 
value of S(kx) in the MS phase, but the vortex posi- 
tions are uncorrelated between channels so there is no 
longitudinal order and the vortex lattice is highly de- 
fected. Since the residual channels have the same period 
as the pinning lattice the transverse scaling in S(k) gives 
long range order. Finally in the ML phase, shown in 
Fig.3(e,f), the vortex positions are disordered and the 
channeling behavior is lost. 

In Fig. 4 we present the central result of this work, 



the dynamic phase diagram as a function of fd and T 
obtained from fixed drive increasing T simulations (de- 
noted by solid circles) and from fixed temperature in- 
creasing drive simulations (denoted by open circles and 
squares). The depinning line is determined from the pro- 
nounced upward curvature in the V(I) curves or the on- 
set of displacements. The MC-MS line is found from 
the saturation of S(ki,) to a minimum value as well as 
the onset of the transverse displacements. The MS-ML 
line is obtained from the point at which the transverse 
S{kx) peak drops to a minimal value and exhibits a scal- 
ing of r\ w 2. The MS-ML line is roughly independent of 
drive above fd = 0.25/o- The MS line shows some curva- 
ture toward lower temperatures for drives near the depin- 
ning line. The MC-MS transition occurs a small amount 
above Tfif = 0.0058 but approaches this value for the 
lower drives. The MS-ML transition line coincides with 
the pinned (zero driving value) melting temperature T^. 
The pinned phase vanishes at T^. 

The onset of these different phases can be understood 
by considering that for higher drives the moving vortices 
spend a great part of their time outside the pinning sites 
because r p < a. For temperatures above the vortices 
enter a molten state while moving between pinning sites 
since they are essentially moving in a clean system. In 
this melted state, the thermal fluctuations overcome the 
vortex-vortex interaction and the correlation of vortices 
in adjacent channels as well as the longitudinal order is 
lost. The vortices start diffusing at random, leading to a 
large increase in d y . Unlike the case of random pinning, 
which induces an additional shaking temperature that 
effectively lowers the temperature at which the vortices 
disorder or melt, vortices moving in a periodic pinning 
array will not experience a shaking temperature. As long 
as T is less than the melting temperature of the non- 
driven clean system T r " p , at any drive the overall moving 
vortex lattice remains ordered. Further there is still a 
pinning effect in the transverse direction. This transverse 
pinning which causes the channeling has been seen in 
simulations with random pinning and is particularly 
large for simulations with periodic pinning arrays |jTo| . 
During the time the vortices are in the pinning sites they 
still feel a transverse pinning force until T > T^, so some 
vortex channeling persists and finite transverse ordering 
appears. Above the pinning is no longer effective so 
channeling and transverse order are both lost. 

A test of the above interpretation is that Tms is related 
to T^f . T^p can be tuned by changing the vortex- vortex 
interaction prefactor A v . As A v is lowered T£f will also 
be lowered. In the inset of Fig. 4 we show the melting 
lines for the clean non-moving system and the moving 
system with pinning (fd = 0.45/o) for T versus A v . As 
A v is increased the T£f and Tms lines increase at the 
same rate. This correlation is encouraging evidence for 
the above interpretation. For the stiffer lattices (high 
A v ) the MC phase is larger. 



In summary, we have studied numerically the melting 
and dynamic phase diagram of a vortex lattice interacting 
with a square periodic pinning array. The melting tem- 
perature for the non-driven pinned system is higher than 
that for the equivalent clean system. For moving vortices 
at low T the vortex lattice moves in correlated ID chan- 
nels and has long range order with some anisotropy be- 
tween the longitudinal and transverse peaks. At a higher 
temperature there is a transition to a Moving Smectic 
phase where longitudinal order is lost while long-range 
transverse order remains. In the Moving Smectic phase 
the channel flow still persists but vortices diffuse between 
channels. At high temperatures the transverse order 
and channeling is also lost. The Moving Crystal-Moving 
Smectic transition corresponds roughly to the melting 
temperature of the zero drive clean system. 

We thank C.J. Olson and R.T. Scalettar for useful dis- 
cussions and critical reading of the manuscript. Funding 
provided by CLC/CULAR. 

-Note Added After submission we became aware of the 
paper by V. Marconi and D. Dominguez fl2|| in which 
they also study the melting of moving vortex lattices in- 
tercting with a periodic substrate in a peroidc Josephson- 
juntion array in which we find several overlapping results. 
They, however, do not find a moving smectic phase which 
is due to their pinning being modeled as an egg-carton 
potential unlike our model in which the space between 
the pinned sites is essential for the smectic phase to oc- 
cur. 
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